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I. INTRODUCTION 



Neutron stars are among the most strongly magnetized objects in the universe fl. Ra- 
dio pulsars are endowed with typical surface magnetic fields of order 10 12 G [2J. A few 
radio pulsars have been found to have significantly higher surface magnetic fields of order 



10 



1.3 



10 14 G [3]. Surface magnetic fields up to 2.4 x 10 15 G have been inferred in soft- 
gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) from both spin-down and 



spectroscopic studies [4j, |5|. Even stronger 
stars, as suggested by various observations 



ields might exist in the interior of these neutron 
6H8( . Duncan and Thompson showed that strong 



magnetic fields up to ~ 10 — 10 G could be generated via dynamo effects in hot newly- 
born neutron stars with initial periods of a few milliseconds [9] leading to the formation of 
strongly magnetized neutron stars thus dubbed magnetars (see e.g. Ref. 10j for a review). 
Numerical simulations confirmed that magnetic fields of order ~ 10 15 — 10 16 G could be 



produced during supernovae explosions due to the magnetorotational instability 



ll|. Huge 



amount of magnetic energy could be occasionally released in crustquakes thus triggering the 
gamma-ray bursts observed in SGRs and AXPs [12]. This scenario has been recently sup- 
ported by the detection of quasiperiodic oscillations (QPOs) in the X-ray flux of giant flares 
from a few SGRs. Some of these QPOs coincide reasonably well with seismic crustal modes 



m 



14j . The huge luminosity variation 



thought to arise from the release of magnetic stresses 
suggests B > 10 15 G at the star surface thus lending support to the magnetar hypothesis [15 ]. 
According to the virial theorem, the upper limit on the neutron-star magnetic fields is of the 
order of 10 18 G [161]. This limit has been confirmed by numerical magnetohydrodynamics 



simulations 



17-1 



In this paper, we study the impact of a strong magnetic field on the equilibrium proper- 
ties of the outer crust of cold non-accreting neutron stars along the lines of Ref. 16j. For 
this purpose we have made use of the most recent experimental atomic mass data comple- 
mented with a theoretical atomic mass table based on the Hartree-Fock-Bogoliubov (HFB) 
method [20|]. In Section [Til we present the microscopic model used to describe the outer 
crust of a magnetar. Results are discussed in Section II III and simple analytical formulae are 
derived in the limit of strongly quantizing fields in Section [IV] 
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II. MICROSCOPIC MODEL OF MAGNETAR CRUSTS 



In the magnetar theory, neutron stars are born with very strong magnetic fields of order 
B ~ 10 16 — 10 17 G which decay on a typical time scale of order > 10 3 years 21] . We assume 



that the magnetic fields are sustained 



ong enough to alter the formation of neutron-star 
crusts. In the model we adopt here [16], the neutron-star crust is assumed to be made of 
"cold catalyzed matter", i.e. matter in its ground state at zero temperature and in a uniform 
magnetic field. The magnetic field mostly affects the outermost region of the crust where 
atoms are supposed to be fully ionised and arranged in a body centered cubic lattice 22]. 
We have determined the equilibrium composition of each layer of the outer crust at a given 
pressure P by minimising the Gibbs free energy per nucleon 

E + P 



9 



n 



(1) 



where £ is the average energy density and n the average nucleon number density. Assuming 
that each layer of the outer crust contains only one nuclear species with proton number Z 
and atomic number A, the average energy density can be expressed as 



£ = n N M'(Z,A) + £ e + £ L 



(2) 



where = n/A is the number density of nuclei, M'(Z,A) their mass (including the rest 
mass of nucleons and Z electrons), £ e the energy density of electrons after subtracting out 
the electron rest mass energy density and El the lattice energy density. The nuclear mass 
M'(Z, A) can be obtained from the atomic mass M(Z, A) after substracting out the binding 
energy of the atomic electrons (see Eq. ( A4) of Ref . [23] ) 

M'{A, Z) = M(A, Z) + 1.44381 x 10 -5 Z 2m + 1.55468 x 10" 12 Z 5 - 35 (3) 

where both masses are expressed in units of MeV. As in Ref. , we will ignore the effects of 
the magnetic field on nuclear masses. Shell correction calculations using the simple Nilsson 
model predict that magnetic fields ~ 10 16 G can change nuclear shell structure hence also 
nuclear masses 



24 



251 ] . However, a very recent study based on fully self-consistent relativistic 



mean-field calculations concludes that significantly hi ghe r fields B <110 17 G are required to 
affect substantially the composition of the outer crust [26]. In Ref. 16|] , the authors used the 
experimental atomic masses from Ref. [2j]] supplemented with the mass model of Ref. 28]. 



3 



In this paper, we have made use of the most recent experimental atomic mass data from a 
preliminary unpublished version of an updated Atomic Mass Evaluation (AME) 2j|. For 
the masses that have not yet been measured, we have employed the microscopic atomic 



mass model HFB-21 of Ref. 
effective nucleon-nucleon interaction 
interaction 



201 ] based on the HFB method using a generalized Skyrme 



30] supplemented with a microscopic contact pairing 



311 ] . The parameters of the Skyrme interaction BSk21 underlying the HFB-21 



model were fitted to the 2149 measured masses of nuclei with N and Z > 8 given in the 
2003 AME 32j. For this it was necessary to add two phenomenological corrections to the 



HFB ground-state energy: (i) a Wigner energy (which contributes significantly only for 
light nuclei or nuclei with N close to Z) and (ii) a correction for the spurious rotational 
and vibrational collective energy. With an rms deviation as low as 0.58 MeV, this atomic 
mass model is well-suited for describing the neutron-rich nuclei found in the outer crust of 
a neutron star. Incidentally, the parameters of the Skyrme interaction were simultaneously 
constrained to reproduce the zero-temperature equation of state of homogeneous neutron 
matter, as determined by many-body calculations with realistic two- and three-nucleon 
forces 33J], from very low densities up to the maximum density found in stable neutron 
stars. For this reason, the Skyrme interaction BSk21 could be reliably extrapolated beyond 
the outer crust thus providing a unified description of all regions of a neutron star. In 
particular, this interaction has been recently used to determine the equation of state of cold 
non-accreting non-magnetized neutron stars_ 34], [35] and has been found to be compatible 



s|3j 



with measurements of neutron-star masses 

In the presence of a strong magnetic field, the electron motion perpendicular to the field 
is quantized into Landau levels (see for instance Ref. [l|). For sufficiently strong fields, 
the electron cyclotron energy becomes comparable to the electron rest-mass energy. This 
happens for B > B c where the critical magnetic field B c is given by 



B r 



1 s 

m e c 
eh 



4.4 x 10" G. 



(4) 

Surface magnetic fields B > B c have been inferred in various kinds of neutron stars jjj, Q]. 
Ignoring electron polarization effects (see e.g. Chapter 4 in Ref. [l| and references therein) 
and treating electrons as a relativistic Fermi gas, the energies of Landau levels (which were 
actually first found by Rabi as early as 1928 37J) are given by 



e v = ^c 2 p 2 z + m\c\\ + 2uB„) 



(5) 



1 

v = n L + - + a 



(6) 



where is any non-negative integer, a = ±1/2 is the spin, p z is the component of the 
momentum along the field, and B+ = B/B c . The electron anomalous magnetic moment 
is small and has been neglected. For a given magnetic field strength B*, The number of 
occupied Landau levels is determined by the electron number density n e 

IB, 



(7) 



u=0 



x e {v) = ^ 1 2 e -\-2vB^ (8) 

where A e = h/m e c is the electron Compton wavelength, 7 e is the electron chemical potential 
in units of the electron rest mass energy, i.e. 

He 



7e 



m P cr 



(9) 



while the degeneracy g u is g v = 1 for v = and g v = 2 for v > 1. 

The electron energy density £ e and corresponding electron pressure P e are given by (see 
e.g. Ref. and references therein) 

x e (v) 



n 9 ^max 

B+m P & 



(2tt)^ 



9u(l + 2uB^ 



u=0 



and 



7 ) O ^max 

B±m e c z 



{2-nfXl 



J29u(l + 2uB^ 



_vTT277a 

x e (u) 



respectively, where 



;Vl + x 2 ± \n(x + VTT 



(10) 



(11) 



(12) 



In the absence of magnetic 
to (see e.g. Chapter 2 in Ref. 

,,2 



lp±(x) = r\ 

iields B = 0, the electron energy density and pressure reduce 
1 



57T 



, ./>(! + 2^) v 7 ! + x 2 r - ln(x r + y/l + x 2 r ) 



and 



8vr 2 A3 



| -s r - 1 ) a/1 + x 2 + ln(x r 



'l + x 



(13) 



(14) 



respectively where x r = fi{?>'K 2 n e ) 1 ^ / \m e c) is the relativity parameter. 

According to the Bohr- van Leeuwen theorem 38j , the lattice energy density is not affected 
by the magnetic field (we neglect here the small contribution due to quantum zero-point 



motion of ions 39] ). For point-like ions arranged in a body centered cubic lattice, the lattice 



energy density is approximately given by 



40] 



£ L = -1.44423Z 2/3 e 2 nf 3 , 



and the associated pressure is 



(15) 



(16) 



The total pressure P is therefore 

P = P e + P L . (17) 

Equations (fTDl) - (fPTj) remain approximately valid at finite temperatures T provided the 
Coulomb coupling parameter r ^> 10 2 and the temperature T is much smaller than the 
electron Fermi temperature Tp. 

III. EQUILIBRIUM COMPOSITION AND EQUATION OF STATE OF MAGNE- 
TAR CRUSTS 



The equilibrium composition of the outer crust of a non-accreting magnetized neutron 
star at T = in a layer characterized by a pressure P is determined by minimizing the 
Gibbs free energy per nucleon 



9 



S + P M'(A,Z) Z 



it A A \ e " 3 n e 

Note that the value of g at equilibrium is simply equal to the neutron chemical potential. 
Starting from the shallowest part of the crust where P ~ 0, we have repeated the calculations 
by increasing the pressure until g equals m n c 2 (m n being the neutron mass) for some pressure 
Pdrip- The present model is not suited for describing the inner regions of the crust because 
neutrons drip out of nuclei for P > Pd r \ v (for a study of the denser regions of strongly 
magnetized neutron stars, see e.g. Ref. 41] for the inner crust and Refs. 42|, |43| for the 
core) . 

Results are summarized in Tables U HlVl For comparison, we have also determined the 
composition of the outer crust in the absence of magnetic field. Note that our results 
shown in Table |V] are slightly different from those given in Table III of Ref. SjJ using the 
same HFB-21 atomic mass model because of our neglect of electron exchange and other 
small corrections that were included in Ref. [34|. Our results for B± = also differ from 



2 4£ L 
m e c H 



(18) 



those obtained previously in Ref. 16] because of the use of more recent experimental and 
theoretical atomic mass data. In particular, the elements 124 Ru and 118 Kr that were found 
in Ref. 16| are now absent, whereas 79 Cu, 80 Ni, 124 Sr and 121 Y are present. Using the latest 
experimental mass tables the outer crust is found to contain 9 nuclides with experimentally 
measured masses vs only 6 in the calculations of Ref. |l6|] . 

For the "weak" magnetic fields prevailing in most pulsars B± < 1, the sequence of equi- 
librium nuclides in their outer crust is the same as that obtained in the absence of magnetic 
fields. However the highest density at which each nuclide can be found is increased, espe- 
cially in the shallow region of the crust where the effects of Landau quantization are the 
most important. For instance, the maximum density at which 56 Fe is found, is raised from 
4.93 x 1CT 9 fm -3 for B* = to 5.60 x 1CT 9 fm -3 for B* — 1. For the strong fields expected 
to exist in magnetars B+ ^> 1, the sequence of equilibrium nuclides is changed. Table |VT] in- 
dicates the magnetic field strength above which a nuclide appears or disappears. Moreover, 
strong magnetic fields tend to prevent neutrons from dripping out of nuclei. The pressure 
at the neutron drip transition thus increases from 4.88 x 10 -4 MeV fm -3 for B* = to 
1.15 x 10 -3 MeV fm -3 for B+ = 2000, as shown in Fig. [TJ Note however that the equilibrium 
nuclide at the neutron-drip point remains 124 Sr in all cases. 

As shown in Fig. [2j the strongly quantizing magnetic field prevailing in magnetar interiors 
is found to have a large impact on the equation of state in the regions where only a few 
Landau levels are filled. In particular, the quantization of electron motion makes the outer- 
most layers of the crust almost incompressible, the density remaining essentially unchanged 
over a wide range of pressures. However, the present model is not well suited for describing 
the surface of the star because of the nonuniformity of the electron gas 44|. In addi 
at finite temperatures thermal effects can considerably change the equation of state 
With increasing density, the effects of the magnetic field become less and less important as 
more and more levels are populated and the equation of state matches smoothly with that 
obtained in the absence of magnetic fields. 



45]. 
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IV. EQUILIBRIUM COMPOSITION AND EQUATION OF STATE OF 
NEUTRON-STAR CRUSTS FOR STRONGLY QUANTIZING MAGNETIC 
FIELDS 



A magnetic field is strongly quantizing if only the lowest level v = is filled. This 
situation occurs whenever the electron number density n e satisfies the inequality (see e.g. 
Chapter 4 in Ref.[l]) 

(19) 



n P < 



V2< 



2 n 3 



where a m = ^hc/eB. Since the average nucleon density is given by n = (A/Z)n e , Eq. ([[9 
can be equivalently expressed as n < tib with 



n B ~ 1.24 x 10- 9 -Bl /2 im.- 3 



(20) 



In the following subsections, electrons will be assumed to fill only the lowest level v = 




200 400 600 800 1000 1200 1400 1600 1800 2000 

B/B 



FIG. 1: Pressure -Pdrip a t the neutron drip transition in the crust of cold non- accreting neutron 
stars as a function of the magnetic field strength = B/B c . The dashed line corresponds to the 
approximate expression (|48p obtained in the strongly quantizing regime for which only the lowest 
level v = is filled. 
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TABLE I: Composition and equation of state of the outer crust of strongly magnetized neutron 
stars with = B/B c = 1. n m i n (n max ) is the minimum (respectively maximum) average baryon 
density (in units of fm -3 ) at which the given nucleus is present. The pressure P ma x (in units of 
MeV fm" 3 ) is the maximum pressure at which the given nucleus can be found. The surface density 
was estimated from Eq. ([271) . 



z 


N 


A 










P 

1 max 




26 


30 


56 


2.50x10- 


10 


5.60x10" 


-9 


4.15x10- 


10 


28 


34 


62 


5.77x10" 


-9 


1.60x10" 


-7 


4.22x10" 


-8 


26 


32 


58 


1.61x10" 


-7 


1.65x10" 


-7 


4.40x10" 


-8 


28 


36 


64 


1.70x10" 


-7 


8.01x10" 


-7 


3.56x10" 


-7 


28 


38 


66 


8.27x10" 


-7 


9.24x10" 


-7 


4.14x10" 


-7 


36 


50 


86 


9.42x10" 


-7 


1.86x10" 


-6 


1.03x10" 


-6 


34 


50 


84 


1.92x10" 


-6 


6.79x10" 


-6 


5.58x10" 


-6 


32 


50 


82 


7.03x10" 


-6 


1.67x10" 


-5 


1.77x10" 


-5 


30 


50 


80 


1.74x10" 


-5 


3.19x10" 


-5 


3.98x10" 


-5 


29 


50 


79 


3.26x10" 


-5 


4.35x10" 


-5 


5.87x10" 


-5 


28 


50 


78 


4.45x10" 


-5 


5.41x10" 


-5 


7.64x10" 


-5 


28 


52 


80 


5.56x10" 


-5 


8.09x10" 


-5 


1.24x10" 


-4 


42 


82 124 


8.37x10" 


-5 


1.22x10" 


-4 


2.07x10" 


-4 


40 


82 122 


1.27x10" 


-4 


1.48x10" 


-4 


2.55x10" 


-4 


39 


82 121 


1.51x10" 


-4 


1.74x10" 


-4 


3.11x10" 


-A 


38 


82 120 


1.78x10" 


-4 


1.95x10" 


-1 


3.53x10" 


-4 


38 


84 122 


1.99x10" 


-4 


2.39x10" 


-4 


4.54x10" 


-4 


38 


86 124 


2.44x10" 


-4 


2.56x10" 


-4 


4.87x10" 


-4 



in all regions of the outer crust. We have found that this assumption is fulfilled whenever 
> 1304. 
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TABLE II: Same as Table Q] for B* = 10. 



z 


N 


A 










P 

1 max 




26 


30 


56 


3.96x10" 


-9 


2.66x10" 


-8 


2.700x10 


-9 


28 


34 


62 


3.06x10" 


-8 


1.85x10" 


-7 


5.67x10" 


-8 


28 


36 


64 


1.93x10" 


-7 


8.14x10" 


-7 


3.75x10" 


-7 


28 


38 


66 


8.40x10" 


-7 


9.14x10" 


-7 


4.36x10" 


-7 


36 


50 


86 


9.31x10" 


-7 


1.85x10" 


-6 


1.05x10" 


-6 


34 


50 


84 


1.91x10" 


-6 


6.74x10" 


-6 


5.59x10" 


-6 


32 


50 


82 


6.99x10" 


-6 


1.67x10" 


-5 


1.77x10" 


-5 


30 


50 


80 


1.73x10" 


-5 


3.19x10" 


-5 


3.98x10" 


-5 


29 


50 


79 


3.26x10" 


-5 


4.36x10" 


-5 


5.88x10" 


-5 


28 


50 


78 


4.46x10" 


-5 


5.42x10" 


-5 


7.64x10" 


-5 


28 


52 


80 


5.56x10" 


-5 


8.00x10" 


-5 


1.24x10" 


-4 


42 


82 


124 


8.36x10" 


-5 


1.23x10" 


-4 


2.07x10" 


-4 


40 


82 


122 


1.26x10" 


-4 


1.48x10" 


-4 


2.56x10" 


-4 


39 


82 


121 


1.50x10" 


-1 


1.74x10" 


-4 


3.12x10" 


-4 


38 


82 


120 


1.77x10" 


-1 


1.95x10" 


-4 


3.53x10" 


-4 


38 


84 


122 


1.98x10" 


-4 


2.40x10" 


-4 


4.55x10" 


-4 


38 


86 


124 


2.44x10" 


-4 


2.57x10" 


-4 


4.88x10" 


-4 



A. Equation of state 



In strongly quantizing magnetic fields, the electron energy density ffTUl) and the electron 
pressure (fill reduce to 

;^+(x e ) - n e m e c , (21) 



respectively, with 



B+rrieC 2 



P<- 



2ir 2 \ 3 e n e 



(22) 



(23) 
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TABLE III: Same as Table U for B± = 100. 



z 


N 


A 










P 

1 max 




26 


30 


56 


6.28x10" 


-8 


2.84x10" 


-7 


2.97x10" 


-8 


28 


34 


62 


2.96x10" 


-7 


1.01x10" 


-6 


5.41x10" 


-7 


28 


36 


64 


1.04x10" 


-6 


1.68x10" 


-6 


1.47x10" 


-6 


36 


50 


86 


1.76x10" 


-6 


2.33x10" 


-6 


2.62x10" 


-6 


34 


50 


84 


2.40x10" 


-6 


7.58x10" 


-6 


7.34x10" 


-6 


32 


50 


82 


7.85x10" 


-6 


1.73x10" 


-5 


1.97x10" 


-5 


30 


50 


80 


1.80x10" 


-5 


3.17x10" 


-5 


4.19x10" 


-5 


29 


50 


79 


3.24x10" 


-5 


4.41x10" 


-5 


6.10x10" 


-5 


28 


50 


78 


4.51x10" 


-5 


5.48x10" 


-5 


7.85x10" 


-5 


28 


52 


80 


5.62x10" 


-5 


8.04x10" 


-5 


1.27x10" 


-4 


42 


82 


124 


8.38x10" 


-5 


1.23x10" 


-4 


2.10x10" 


-4 


40 


82 


122 


1.27x10" 


-4 


1.49x10" 


-4 


2.58x10" 


-4 


39 


82 


121 


1.51x10" 


-4 


1.75x10" 


-4 


3.14x10" 


-4 


38 


82 


120 


1.78x10" 


-4 


1.95x10" 


-1 


3.55x10" 


-1 


38 


84 


122 


1.99x10" 


-4 


2.40x10" 


-4 


4.57x10" 


-4 


38 


86 


124 


2.44x10" 


-4 


2.57x10" 


-1 


4.90x10" 


-4 



The electron chemical potential can be obtained from 



7e 



xf 



(24) 



In the upper layers of the crust where i e C 1, the electron pressure (l22l) is approximately 
given by 



P e rti -m e c 2 nl 
3 



2tt 2 A 3 
5. 



(25) 



Substituting Eq. fl25|) into Eq. f|T7|) with P = using Eq. f fT6|) yields the average density at 
the surface of a cold non-accreting magnetar 16j 

2 n 3/5 

(26) 



A, 



1.44423Z s 2/3 e 2 



P, 



2tt 2 A 3 



with Z s and A s the proton number and the charge number of the equilibrium nuclide at the 
surface. Considering that the surface of a neutron star is made of iron with Z s = 26 and 
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TABLE IV: Same as Table U for = 1000. 



z 


N 


A 


1T"min 








P 

1 max 




26 


30 


56 


9.96x10" 


-7 


2.62x10" 


-6 


1.98x10" 


-7 


28 


34 


62 


2.71x10" 


-6 


1.10x10" 


-5 


6.23x10" 


-6 


28 


36 


64 


1.14x10" 


-5 


1.42x10" 


-5 


1.01x10" 


-5 


38 


50 


88 


1.45x10" 


-5 


1.55x10" 


-5 


1.16x10" 


-5 


36 


50 


86 


1.60x10" 


-5 


2.60x10" 


-5 


3.21x10" 


-5 


34 


50 


84 


2.69x10" 


-5 


3.88x10" 


-5 


6.81x10" 


-5 


32 


50 


82 


4.03x10" 


-5 


5.22x10" 


-5 


1.16x10" 


-1 


30 


50 


80 


5.43x10" 


-5 


6.54x10" 


-5 


1.69x10" 


-4 


29 


50 


79 


6.68x10" 


-5 


7.32x10" 


-5 


2.03x10" 


-4 


28 


50 


78 


7.48x10" 


-5 


7.92x10" 


-5 


2.28x10" 


-4 


28 


52 


80 


8.12x10" 


-5 


9.03x10" 


-5 


2.83x10" 


-4 


42 


82 


124 


9.37x10" 


-5 


1.07x10" 


-4 


3.70x10" 


-4 


40 


82 


122 


1.10x10" 


-4 


1.16x10" 


-1 


4.09x10" 


-4 


39 


82 


121 


1.18x10" 


-4 


1.79x10" 


-4 


4.59x10" 


-1 


38 


82 


120 


1.82x10" 


-4 


2.17x10" 


-1 


5.04x10" 


-4 


38 


84 


122 


2.21x10" 


-4 


2.74x10" 


-4 


6.20x10" 


-4 


38 


86 


124 


2.78x10" 


-4 


2.92x10" 


-4 


6.57x10" 


-4 



A s = 56 leads to 

n s ~ 2.5 x KT 10 ^/ 5 frrT 3 . (27) 

This simple formula shows that the stronger the magnetic field is, the higher is the surface 
density. It should be stressed however that Eq. (1271) provides only an approximate estimate 
of the surface density because for sufficiently strong fields the condition i e C 1 is not 
fulfilled. Moreover, the present model is not strictly valid at the surface of a neutron star, 
as mentioned earlier. Using Eq. ( 1231) and ( )26l) . the condition i e < 1 at the neutron star 
surface translates to 

2rr 2 / he \ 3 , 

< — -) ~2.5xl0 4 . 28 

Z 2 S Vl-44423e 2 y V ; 

In the dense region of the outer crust where n 3> n s , the lattice pressure (TLB"]) is negligible 
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TABLE V: Same as Table U for = 0. 



z 


N 


A 










P 

1 max 




26 


30 


56 







4.93x10" 


-9 


3.36x10- 


10 


28 


34 


62 


5.09x10" 


-9 


1.59x10" 


-7 


4.20x10" 


-8 


26 


32 


58 


1.60x10" 


-7 


1.65x10" 


-7 


4.39x10" 


-8 


28 


36 


64 


1.70x10" 


-7 


7.99x10" 


-7 


3.55x10" 


-7 


28 


38 


66 


8.26x10" 


-7 


9.22x10" 


-7 


4.13x10" 


-7 


36 


50 


86 


9.42x10" 


-7 


1.86x10" 


-6 


1.03x10" 


-6 


34 


50 


84 


1.92x10" 


-6 


6.79x10" 


-6 


5.57x10" 


-6 


32 


50 


82 


7.05x10" 


-6 


1.67x10" 


-5 


1.77x10" 


-5 


30 


50 


80 


1.74x10" 


-5 


3.18x10" 


-5 


3.98x10" 


-5 


29 


50 


79 


3.26x10" 


-5 


4.35x10" 


-5 


5.87x10" 


-5 


28 


50 


78 


4.46x10" 


-5 


5.42x10" 


-5 


7.64x10" 


-5 


28 


52 


80 


5.57x10" 


-5 


7.99x10" 


-5 


1.24x10" 


-4 


42 


82 


124 


8.36x10" 


-5 


1.23x10" 


-1 


2.07x10" 


-4 


40 


82 122 1.27x10" 


-4 


1.48x10" 


-1 


2.55x10" 


-4 


39 


82 121 


1.51x10" 


-4 


1.74x10" 


-1 


3.11x10" 


-4 


38 


82 120 1.78x10" 


-4 


1.95x10" 


-4 


3.53x10" 


-4 


38 


84 122 1.99x10" 


-4 


2.39x10" 


-4 


4.54x10" 


-4 


38 


86 124 2.44x10" 


-4 


2.56x10" 


-4 


4.86x10" 


-4 



and the total pressure is approximately given by 



9 9 7r 2 A^ 
P ~ P e K, m e c n e - 



Inverting this equation yields 



A ( PB+ 
n = — 



1/2 



(29) 



(30) 



Z \m e c 2 7r 2 Ag 

Interpolating between the shallow and the deep regions of the outer crust, the density n 
in a layer at pressure P can be approximately expressed as 

T 



nsinJ 1 + 



(31) 



13 



TABLE VI: Magnetic field strength B± = B/B c for the appearance (+) or the disappearance (-) 
of a nuclide in the the outer crust of a cold non-accreting neutron star. 



Nuclide 


B, 


58 Fe(-) 


9 


66 M (.) 


67 


88 Sr(+) 


859 


126 Ru(+) 


1118 


128 Pd(+) 1120 


78 Ni(-) 


1120 


80 Ni(-) 


1250 


64 Ni(-) 


1668 


79 Cu(-) 


1791 


130 Cd(+) 1804 


132 Sn(+) 


1987 



where 



2 n]^\l (Z 



1.82 x 1Q- 11 Bl /5 { j ) MeVfm" 3 



(32) 



As illustrated in Fig. [31 the analytical representation fl3T|) yields a fairly good fit to the 
equation of state obtained from the full minimization of the Gibbs free energy g. The 
typical error is found to be less than 11% for B± > 10 in any region of the outer crust where 
the condition (flUI) holds. 

At densities n > % or equivalently at pressures P 3> Pb, many Landau levels are 
populated so that the quantization effects disappear and the properties of the crust are 
almost unaffected by the magnetic field. According to Eqs. (1201) and (!30|) . the pressure Pb 
is given by 

Pb = $^. (33) 



2tt 2 \l 



B. Composition 

The results about the crustal composition presented in Section [III] can be qualitatively 
understood using a simplified atomic mass formula. Neglecting Coulomb and surface con- 
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tributions, the mass of a nucleus with Z protons and A nucleons is given by 
M'(A, Z) = A(a v + J(l - 2y e f + m u c 2 ) + Zm e c 2 , 



(34) 



where y e = Z/A is the electron fraction, a v is the binding energy of symmetric nuclear 
matter, J the symmetry energy and m u the atomic mass unit (ignoring here the small 
difference between neutron and proton masses). For the HFB-21 nuclear mass model that 



we consider here 



20|, a v 



-16.053 MeV and J = 30 MeV. Dropping the lattice energy 



density El (which is a small correction to the total energy density E), the Gibbs free energy 
per nucleon as given by Eq. (fl8|) reduces to 

g = a v + J(l - 2y e f + m u c 2 + y e /i e . (35) 

Minimizing Eq. (135]) for a given pressure P ~ P e (i.e. /i e fixed) and treating y e as a continuous 
variable yield 



Ve 



1 _ /fe 

2 8J 



(36) 




n [fm ] 



FIG. 2: Pressure P vs average nucleon number density n in the outer crust of a cold non-accreting 
neutron star for different magnetic field strengths B+. The filled squares indicate the points above 
which the lowest level v = is fully occupied. 
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Using Eqs. and (129"]) leads to 



se= 5 i_ v£ 



(37) 



where 



P, 



4£LJ 2 



neu 2 \ S 2 



(38) 



For comparison, using Eqs. ( !T3|) and (fl~4"|) the electron fraction in the absence of magnetic 
field is approximately given by 

l/4> 



where 



""=2 l 1 "^ 



4 J \ 4 m e c 2 



m e c 2 J 127r 2 Ag 



(39) 



(40) 



10" 



10" 



10" 



10" u 

> 

io" 7 



10" 



10" 



10 



10 




10 



-6 



analytical representation 
full minimization of g 



10 

-3 

n [fm J 



10" 



FIG. 3: Pressure P vs average nucleon number density n in the outer crust of a cold non-accreting 
neutron star for = 1400 from the full minimization of the Gibbs free energy (solid line) and 
from the analytical representation (|3ip (dashed line). 
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assuming P ~ P e and x r ^> 1. Introducing the isospin asymmetry parameters 77 = 1 — 2y e 
and 77° = 1 — 2y®, their ratio is given by 



P 



1/4 / p \ 1/4 
neu 



(41) 



T7 U \ P / VP 

'/ \ J ncu / \ J ncu, 

As will be shown in the next section P <C P neu in any region of the outer crust. Noting that 

JO \ 1/4 



pu 



J 



< 2 



(42) 



Pneu^ V m e C 2 (3P,) V4 

for P + > 1304 (strongly quantizing field), we find that rj < rf. In other words, nuclei in the 
outer crust of a magnetar are more symmetric than those found in the outer crust of a weakly 
magnetized neutron star at the same pressure. This conclusion is confirmed by numerical 
calculations using the experimental and HFB-21 atomic masses, as shown in Fig. HJ 




3 4 5 
P[10" 4 MeVfm 3 ] 



FIG. 4: Electron fraction y e vs pressure P in the outer crust of a cold non-accreting neutron star for 
B* = 1400 (solid line) and for B* = (dashed line). Note that in both cases the equilibrium nuclide 
at the bottom of the outer crust is 124 Sr even though the neutron drip pressures are different. 
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C. Neutron drip transition 



With increasing pressure, the crustal matter becomes more and more neutron rich, as 
shown by Eq. (j57|) . The pressure P neu can be interpreted as the pressure at which nuclei 
would convert to neutron drops. In reality, neutrons start to drip out of nuclei above some 
pressure Pdrip which is lower than P neu - This transition occurs when the neutron chemical 
potential (i n exceeds the neutron rest mass energy. As can be seen in Tables HHV] the 
equilibrium nucleus at the neutron drip point is independent of the magnetic field strength 
and is found to be 124 Sr for the HFB-21 atomic mass model considered here. The reason is 
the following. Equilibrium with respect to weak interaction processses requires 

Hp + He = Hn , (43) 

where \i v is the proton chemical potential. However at equilibrium, the neutron chemical po- 
tential coincides with the Gibbs free energy per nucleon g. Neglecting the small contribution 
of the lattice energy density, Eq. ffT8]) leads to 

2 A ( M'(A,Z)\ 
H e « m e c 2 + -U, A^J ■ (44) 

Substituting the neutron drip value of the neutron chemical potential fi n = m n c 2 in Eqs.(j43l) 
and (jH]), we find 

2 ^ A fM'(A,Z) 



pLp - m p c = Qn :l3 + — ^ m n c J , (45) 

where m p is the proton mass and Q n ,/3 is the beta decay energy of the neutron. The quantity 
on the left-hand side of Eq. (1451) is approximately equal to the opposite of the one-proton 
separation energy. The equilibrium nucleus at neutron drip is therefore uniquely determined 
by atomic masses. Using the two-parameter mass formula discussed in Section HVBt we find 
that the proton fraction is approximately given by 



^-2V 1 + T- <«> 

Substituting the values of a v and J from the atomic mass model HFB-21 in Eq. (T46l) yields 
a fairly good estimate of the proton fraction of 124 Sr with an error of about 11% only. Note 
however that without Coulomb and surface terms in the mass formula, it is not possible to 
determine Zdrip and A& V ip separately. 
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Eq. (jUJ) shows that the electron chemical potential at neutron drip is independent of the 
magnetic field and is given by 



p = m eC 2 + (m n c 2 - M/(A f p,Zdrip) ) . (47) 

•^drip \ ^drip / 

For 124 Sr, we find fif ip ~ 26 MeV. Using Eqs. (EH}, (E3}, (ESJ and (ED implies that the 
pressure at the neutron drip point increases linearly with the magnetic field strength (in the 
strongly quantizing regime) as shown in Fig. [I] and is given by 

_ m e c 2 ( 7 drip ) 2 , , 

P dri P - A 3 ^ B*. (48) 

The corresponding baryon density is given by 

A, . ,-ydrip 

— itM B - (49) 

Using the two-parameter mass formula yields 



P.lnp = Pn,u ( 1 - \/l + 7 ) < Pn,u- (™) 



— 1 



^V 1 + 7 O-M)- (51) 

Electrons fill only the lowest level v = in any region of the outer crust provided Pb > 
-Pdrip- Using Eqs. (133]) and fj48l . we find that this condition is equivalent to B± > B+ Tip with 

5f P = ^(7 rinp ) 2 - (52) 

This estimate could have been immediately obtained from Eq. (|Sj) requiring x\ > 0. For 
124 Sr, we find 5f ip ~ 1300. 

In the absence of magnetic fields, the neutron drip pressure and baryon density (in the 
ultrarelativistic regime x r ^> 1) are approximately given by 

^ ~ A3 12tt 2 (53) 

o ^ A dvip ( 7e drip ) 3 
drip ~ Z drip 3vr 2 A3 ' 
respectively. Using Eqs. (HE}, (152"} and ( 155|) leads to 

-^"p _ 3 B * > 3 

-'drip z -D* z 

This shows that the neutron drip transition occurs at a higher pressure in a magnetar than 
in a weakly magnetized neutron star. 
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(54) 



D. Elastic properties 



Because a sufficiently strong magnetic field changes the composition of the outer crust of 
a neutron star, it can also have an impact on the crustal properties. In view of the recent 
detection of QPOs in the X-ray flux of giant flares from SGRs, a particularly important 
property of strongly magnetized neutron star crusts is the shear modulus which determines 
the frequencies of torsional oscillations. 

We have calculated the "effective" shear modulus S of the outer crust, assuming that it 



is made of a body-centered-cubic lattice polycrystal, using the following expression 



46| : 



Z 2 e 2 

S = 0.1194ri N — - , (56) 



where i?N is the ion-sphere radius defined by 



3 ^ 1/3 



As shown in Fig. El the effective shear modulus of the outer crust of a neutron star can be 
enhanced by the presence of a strong magnetic field. 

E. Global structure 

In the outer crust of a non-rotating neutron star of gravitational mass Ai and circumfer- 
ential radius R, the general relativistic equations can be approximately written as (see e.g. 



m 

dp 

where g s is the surface gravity defined by 



9sp, (58) 



GM a r„\-i/2 
9. = ^r(i-4) (59) 



R 2 V R, 

z is the depth below the surface, r g = 2GM/c 2 is the Schwarzschild radius and p ~ nm u 
is the mass density. The gravitational mass AA4 M. contained in the outer crust is 
approximately given by 

AM » g - l) . (60) 

Comparing Eqs. ( 148]) and (|53|) shows that the crustal mass for neutron stars endowed with 
strongly quantizing magnetic fields is larger than that of weakly magnetized neutron stars 
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with the same mass and radius, and is given by 



AM = 2fi 3 m e ( 7e drip ) 2 i?, fR _ 1 
7r XI \r g 



(61) 



Likewise the magnetic field increases the baryonic mass contained in the outer crust, which 
is approximately given by 

(62) 



AM, 



B 



1 - -f AM . 
R 



On the contrary, the depth z below the surface where neutron drip occurs and which therefore 
delimits the boundary between the outer and inner crusts, does not depend on the magnetic 
field strength. Indeed, in the absence of magnetic fields assuming that the main contribution 
to the pressure is due to ultra-relativistic electrons, i.e., P oc p 4//3 , Eq. ( 158|) can be easily 



solved, leading to 



24 1 



4drip C 



M&pfl lR fR 



- l 



(63) 








3 4 5 
P[10" 4 MeVfm 3 ] 



FIG. 5: Effective shear modulus S vs pressure P in the outer crust of a cold non- accreting neutron 
star for = 1400 (solid line) and for = (dashed line). Note that the neutron drip pressures 
are different in the two cases. 
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which can be expressed as 



n 2m e 7f ip ?/f i PR \R f R 



m v , \i r g \r g 



z°= e '° °° kE- ^-1 , (64) 



where we have used Eqs. (153)) and (151)1 . In the presence of a strongly quantizing magnetic 
field, Eq. (129)) shows that the pressure varies approximately as P oc p 2 . Solving Eq. (l58)) thus 
yields 

AP dlip R [R (R 1( (65) 



Using Eqs. flU} and (gS} leads to Eq. (JM)) so that z = z°. 
V. CONCLUSION 

We have calculated the composition and the equation of state of the outer crust 
of cold non-accreting neutron stars endowed with very strong magnetic fields of order 
B ^> m 2 e c'/{eh) ~ 4.4 x 10 13 G, as measured in soft-gamma repeaters, anomalous X-ray 
pulsars and even in a few radio pulsars (3), 5^ For this purpose, we have made use of the 



most recent experimental atomic mass data 



Bogoliubov atomic mass model 



20] 



29j complemented with the latest Hartree-Fock- 



The Landau quantization of electron motion due to the strong magnetic field is found 
to have a significant impact on the neutron-star crust properties: i) it changes the crustal 
composition (the sequence of equilibrium nuclides being different than that found in weakly- 
magnetized crusts, as summarized in Table IVip and iii) it makes the matter less neutron rich 
as shown in Fig. H] and tends to prevent neutrons from dripping out of nuclei (the pressure 
at neutron drip increasing with B as shown in Fig. [TJ. As a consequence the presence of a 
strong magnetic field can have an impact on the crustal properties like the shear modulus, as 
shown in Fig. |5j These results may have implications for the interpretation of quasiperiodic 
oscillations observed in soft gamma-ray repeaters. Likewise, other crustal properties such as 
the thermal and electric conductivities could be affected. The present results might therefore 
also impact the thermal and magnetic field evolution of magnetars. This warrants further 
studies. 

The outer crust of a magnetar is also found to be much more massive than the outer crust 
of a weakly magnetized neutron star with the same gravitational mass M. and circumferen- 
tial radius R. This implies that the contribution of magnetars to the galactic enrichment in 
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nuclides heavier than iron from the rapid neutron capture process (r-process) of nucleosyn- 
thesis following the ejection and the decompression of crustal material 47| . could be much 
more important than previously thought. 
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